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Abstract. A new set of evolutionary synthesis spectra are presented for Simple Stellar Populations (SSPs) 
covering ranges in metallicity from 0.02 < Z/Zq < 2.5 and ages from 4 ■ 10 6 yr < t < 16 Gyr. They are based 
on the most recent isochrones from the Padova group that extend earlier models by the inclusion of the thermal 
pulsing AGB phase for stars in the mass range 2 Mq < m < 7 Mq in accordance with the fuel consumption 
theorem. We show that with respect to earlier models, inclusion of the TP- AGB phase leads to significant changes 
in the (V — I) and (V — K) colors of SSPs in the age range from 10 8 to > 10 9 yr. Using model atmosphere 



spectra from Lejeune et al. (1997, 1998), we calculate the spectral evolution of single burst populations of various 
metallicities covering the wavelength range from 90 A through 160 /mi. Isochrone spectra are convolved with 
filter response functions to describe the time evolution of luminosities and colors in Johnson, Thuan & Gunn, 
Koo, HST, Washington and Stromgren filters. 

The models and their results are not only intended for use in the interpretation of star clusters but also for 
combination with any kind of dynamical galaxy formation and/or evolution model that contains a star formation 
criterion. Moreover, the evolution of these single burst single metallicity stellar populations is readily folded with 
any kind of star formation - and eventually chemical enrichment - history to describe the evolutionary spectral 
synthesis of composite stellar populations like galaxies of any type with continuous or discontinuous star formation. 
For these latter purposes we also present the time evolution of ejection rates for gas and metals for two different 
Initial Mass Functions (IMFs) as well as cosmological and evolutionary corrections for all the filters as a function 
of redshift for < z < 5 and two different cosmologies. 

Extensive data files are provided in the electronic version, at CDS, and at our above www-address. 
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1. Introduction 

SSPs are single burst single metallicity stellar populations 
like star clusters with all their stars having the same age 
and the same metallicity. Evolutionary synthesis models 
describe the evolution of these SSPs in terms of photom- 
etry and spectra. They are of interest for several reasons. 
1. They allow to test a set of stellar evolutionary tracks 
against observations of star clusters from young open clus- 
ters to old Globular Clusters (GCs). 2. They are used for 
the interpretation of star clusters formed or forming in 
strong starbursts like those accompanying galaxy inter- 
actions or mergers. Being formed from the pre-enriched 
ISM of the merging galaxies, these young star clusters 
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may cover a wide range in metallicity. 3. Evolutionary 
models for individual stellar generations can directly be 
folded - in terms of luminosity or spectral evolution - 
with any Star Formation History (SFH) to describe the 
composite stellar population of a galaxy. Or, the other 
way round, a complete library of spectra for SSPs of var- 
ious ages and metallicities provides a theoretical basis for 
population synthesis models. Population Synthesis models 
search for a particular linear combination of input spectra 
that gives best agreement with the observed spectrum of a 
galaxy with unknown SFH. 4. Evolutionary synthesis re- 
sults for SSPs can be used in any kind of dynamical galaxy 
formation and/or evolution code that involves a star for- 
mation criterion. An example is presented in Contardo et 
al. ( 1999 ), where - in the context of a cosmological struc- 
ture formation scenario - an individual galaxy is followed 
from the onset of Star Formation (SF) in a couple of 
subgalactic fragments to a state where it resembles the 
Milky Way in many respects. Spatially resolved observ- 
able quantities are obtained from a superposition of SSP 
spectrophotometric evolutionary model results. 
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A basic difficulty for single burst models as opposed to 
models for galaxies - which always have extended and con- 
tinuous SFHs - is the discreteness of the stellar mass spec- 
trum for which evolutionary tracks have been calculated. 
Lifetime differences between stars for which stellar evolu- 
tionary tracks are available are of the order of > 5 • 10 8 
yr, i.e. much longer than the age spread among the stars 
in SSP models with burst durations of < 10 5 yr and also 
much longer than a typical timestep in the evolutionary 
synthesis model. In late evolutionary phases, this causes 
large numbers of stars to collectively enter and leave dis- 
tinct phases of evolution with the effect of significant os- 
cillations in colors and magnitudes. These oscillations in 
the photometric properties have been cured in different 
ways. E.g., they can be smoothed a posteriori. This has 
the advantage that the intrinsic uncertainty of the method 
remains clearly visible. A smooth spectral evolution, how- 
ever, is difficult to obtain. A second possibility to start 
from a discrete set of stellar tracks and obtain a reason- 
ably smooth photometric evolution is the use of a Mont e 
Carlo method as described e.g. by Kurth et al. ( 1999| ) . 
This method works very well for a photometric descrip- 
tion but requires a huge number of Monte Carlo stars and 
timcstcps to yield a smooth spectral evolution. The third 
possibility is the use of isochrones obtained on the basis 
of a grid of stellar evolutionary tracks by interpolating 
between equivalent evolutionary stages. Once the defini- 
tion of equivalent evolutionary stages is taken for granted, 
this method provides a smooth evolution of SSPs, both in 
terms of photometry and spectra. 

Here we present a new set of evolutionary synthesis 
models for the photometric and spectral evolution of SSPs 
of various metallicities 0.02 < Z/Zq < 2.5 over a Hubble 
time, based on the latest Padova isochrones (cf. Sect. 2.2). 
In this se nse, o ur models are complementary to Leitherer 
et al. 's ( |1999D Starburst99 models which provide a very 
detailed description of the earliest evolutionary phases 
(< 1 Gyr) of instantaneous burst populations (— SSPs). 
In Sect. H we present our models and the input physics we 
use. Sect. | gives the spectral evolution of SSPs for 5 dif- 
ferent metallicities from 10 6 yrs to 10 10 yr and a discussion 
of the effects of metallicity on the time evolution of the 
spectra. In Sect. |], we present the evolution of photomet- 
ric quantities for a number of widely used filter systems. In 
Sect. §, we compare our results to earlier models and GC 
data. In Sect. ^, we present cosmological and evolutionary 
corrections for our SSP models in 2 different cosmologies 
in a way that can immediately be used by observers or in- 
cluded into dynamical galaxy formation and evolutionary 
codes. 

Only a few selected figures are presented here, the full 
tables for the time evolution of spectra, luminosities, col- 
ors, ejection rates, etc. of SSPs of all metallicities as well 
as the corresponding evolutionary and cosmological cor- 
rections as a function of rcdshift can be found at our web- 



SFHs, as well. In this paper, however, we only present SSP 
model results. 



2. Models and Input Physics 

2.1. Input physics: IMF 

Results in this paper are presented for two different forms 
of the stellar IMF. While a Salpeter IMF is parametrised 



as usual <i>(m) 



-2.35 



over the range of stellar masses 



m from mi = 0.15 M to m u w 85 M Q , a Scalo IMF 
is used in the form <E>(m) ~ m~ x with x = —1.25 for 
m < 1 M Q , x = -2.35 for 1 Mq < m < 2 M , and 
x = —3.00 for to > 2 Mq. Normalization in both cases is 
/mi" ®( m ) ' m ' dm = 1. This implies that - as compared 
to a Salpeter IMF - the Scalo IMF has more stars in the 
mass range 0.5 Mq < m < 6.5 Mq and less stars below 
0.5 Mq and above 6.5 M . 

2.2. Input physics: Padova isochrones with TP-AGB 

The results from our evolutionary synthesis model are 
based on isochrones from November, 18th, 1999 of the 
Padova group (for a description and discussion of pa- 
rameters for these isochrones with moderate core over- 



shooting sec Bcrtclli et al. 1994 and Girardi et al. 2000 ). 
Comparison of these isochrones with CMDs of resolved 
star clusters in the Milky Way and LMC were e.g. pre- 



page ( ht t p : / / w w w . uni-s w . g wdg . do / ^ gale v / ) 



Our isochronc synthesis code is able to describe the 
evolution of stellar populations in galaxies with extended 



sented by Chiosi et al. ( |1992|) and Vallenari et al. ( pL994j ) . 
We used the following data from their files: mass, log(TE), 
A/boI, arid My in their time evolution. We do not use the 
colors and other luminosities also contained in their files, 
because we obtain them from the spectra we calculate. 

For the use with single short burst stellar populations, 
isochrones have a basic advantage over the stellar evolu- 
tionary tracks they are based on. Because tracks are only 
available for typically 50 (or less) stellar masses, drastic 
changes in the luminosity evolution of an SSP occur when 
all stars of a certain mass collectively climb on the red gi- 
ant branch or die. Between reasonably defined equivalent 
evolutionary stages of two stellar masses with evolution- 
ary tracks available, isochrones are interpolated in terms 
of 

logL = log (L (mi, /j, r^)) 
logT eff = log (T off (TOi, /i, t m )) 

where m\ is the initial mass for track i, /i the defined phys- 
ical phase of evolution and r M the relative duration of this 
phase. 

It is then possible to interpolate both in luminosity 
and temperature between two adjacent masses. For the 
interpolated masses the magnitudes and colors are then 
obtained by translating the theoretical luminosities and 
effective temperatures with tables of bolometric correc- 
tions and colors. 

We supplemented these Padova isochrones by low mass 
stars in the mass range from 0.15 — 0.4 5 Mq from the 
calculations of Chabrier & Baraffe ( 1997 ) that include a 
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new description of the interior of low mass objects and 
use non-grey atmospheres (see Kurth et al. 199£ for de- 
tails). The Padova group recently added low mass-stars to 
their isochrones themselves (Girardi et al. 200C) but did 
not compare their low mass stars to those of Chabrier 
& Baraffe. The cumulative contribution of stars below 
0.5 Mq, however, to the integrated light is very small, 
anyway. 
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Fig. 1. Time evolution of (V — K) and (V — I) colors 
(as a function of logarithmic age [yr]) for SSPs based on 
isochrones, i.e. including the TP-AGB phase (solid 
lines), and based on stellar evolutionary tracks that did 
not yet include the TP-AGB phase (dashed lines), both 
for solar metallicity. 



An important feature of the 1999 Padova isochrones 
is that they include a semi-analytical description of the 
thermal pulsing asymptotic gia nt bra nch (TP-AGB) phase 
(cf. Groenewegen & de Jong 1993 ), which was not yet 
inclu ded in the ear lier Padova tracks (c f. Bressan et al. 
1993| , Fagotto et al. |l994a| , |l994b| , |l994d) . 



TP-AGB stars are rare species in single burst popu- 
lations, they are, however, very luminous and very cool 
(Lancon 1999 ). At ages between 0.1 and 1 Gyr, i.e. when 
the red supergiants have died already and the red giants 
are not important yet, TP-AGB stars account for 25 to 40 



% of the bolometric light of an SSP and for 4 to 6 % of 
the light emitted in the K— band (cf. Chariot 199£). 



Inclusion of the TP-AGB phase has strong effects on 
some of the colors. It makes the broad band colors redder 
in the age range from 10 s to a few 10 9 yr. This effect is 
strongest in (V — K) and can amount to more than 1 mag 
in (V — K) at solar metallicity as seen in Fig. [j], where 
we plot the (V — K) color evolution for an SSP based 
on isochrones including the TP-AGB in comparison with 
an SSP from Kurth et al. ( 1999 ) based on the Padova 
stellar evolutionary tracks that did not include the TP- 
AGB phase. The effect is weaker at lower metallicities and 
for colors at shorter wavelengths. 

The effect of the TP-AGB phase has been discussed 



by the Padova group before (Girardi et al. 1998) and the 
good agreement with J— and K— CMDs of resolved young 



open clusters is shown in Vallenari et al. (1999). Rich et 
al. (2001) find excellent agreement of the recent Padova 



isochrones with CMDs in UBV and / for the LMC cluster 
NGC 2121 which, at [Fe/H] = -0.6 and an age of 3.2 Gyr, 
is a very good test case for the importance of the TP-AGB 
phase. 

The effect of the TP-AGB phase is also significant in 
(V — I) at metallicities > 0.5 Zq and has important im- 
plications for the age determination on the basis of colors 
observed for bright star clusters in interacting and/or star- 
burst galaxies. E.g. in our analysis of the sample of young 
star clusters in the merger remnant NGC7252 which show 
a mean (V — I) = 0.83 ± 0.25 we derived a mean age 
of 1.37 Gyr assuming a metallicity ~ 0.5 Zq using 
SSP models based on Geneva stellar evolutionary tracks 
(Fritze - v. Alvensleben & Burkert 1995) very simi lar to 



what is found with the models of Kurth et al. ( 1999 ) who 
use Padova tracks. Including the TP-AGB phase as in the 
Padova isochrones used here, however, results in a mean 
age of 8.7 • 10 8 yr, i.e. only 64% of the age obtained with- 
out the TP-AGB phase. Most young star cluster systems 
observed with HST, in fact, have ages in the range of 10 s 
to 10 9 yr and metallicities > 0.5 Zq where the effect of 
the TP-AGB phase on (V — /) is strong. Hence, all their 
age determinations on the basis of (V — I) in the literature 
are subject to this effect. 



2.3. Input physics: spectra 

We use the library of model atmosphere spectra from 
Lejeune et al. (1997, 1998) which is based on model at- 
mospheres from Kurucz (see e.g. Kurucz 1992), Fluks et 
al. ( |1994| ), and Bessell et al. fll989| , |1991| ). This library 
comprises the full range of metallicities covered by the 
isochrones, is very complete in terms of stellar effective 
temperatures from T ff = 2800 — 47500 K, and gravi- 
ties, and has a long wavelength coverage from 90 A to 160 
/im with good spectral resolution. Its spectra have partly 
been corrected to provide good agreement with the ob- 
served colors of all types of stars from U through K (cf. 
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Lejeune et al. ). Stars with an effective temperature above 
50000 K are described by pure black body spectra. 

All the coolest AGB stars are long period variables 
which have different spectral type - temperature - color 
relations and, in particular, strongly increased molecular 
absorption features (e.g. H2O around 1.4 and 1.9 /im, VO 
at 1.05 /im, . . . ) as compared to non-pulsating stars 
(Langon 1999). Lanqon & Wood (1997) compiled an IR - 
library of cool star spectra that was combined with the 
one of Lejeune et al. (1998) that we use here and im- 
plemented into models by Fioc & Rocca - Volmerange 
(1997). These models show that the inclusion of observed 



AGB spectra is very important for a description of specific 
spectral features in the near and mid-IR of intermediate 



age SSPs. Bressan et al. (199S) investigate in detail the 
effects of dust in the envelopes of AGB stars on the inte- 
grated IR properties of intermediate age SSPs. They con- 
sider both the absorption of light in the expanding dust 
shells of mass-losing stars at optical wavelengths and its 
thermal reemission in the IR. Differences with respect to a 
standard SSP without these dust absorption and emission 
features appear beyond a few /jm where the contribution 
of the brightest AGB stars to the integrated light of the 
stellar population is large. Bressan et al. find a decrease 
of the flux between 1 and 3 /im because the brightest 
stars are heavily obscured by their circumstellar envelopes 
and thermal dust emission at longer wavelengths > 4 /im, 
reaching an order of magnitude at 10 /im over the pure 
photospheric models. 

The UV, optical and NIR regions we are studying in 
this paper are not affected by the expanding dust shells 
surrounding TP-AGB stars. We caution, however, that 
our models do not include the above-mentioned empiri- 
cal AGB spectra yet, nor the absorption and reemission 
of light by their dusty envelopes, and, hence, should not 
be used for the interpretation of fluxes or spectral features 
of intermediate age star clusters beyond the K— band. In 
stellar populations with extended SFHs the AGB star fea- 
tures discussed above (molecular absorption bands, ex- 
tinction and reemission by surrounding dust shells) are 
strongly overwhelmed by the emission of coexisting red 
supergiant and giant stars. 

2.4. Isochrone synthesis model 

Our isochrone synthesis code for SSPs starts from a gas 
cloud of given mass and metallicity and then forms a pop- 
ulation of stars with a given SFR during the first timestep 
of the calculation. It calculates a variety of observables at 
any evolutionary timestep, typically 4 • 10 6 yr: 

— masses in gas and stars 

— ejection rates of gas and metals by dying stars 

— spectra from 90 A to 160 /im 

— luminosities, mass-to-light - ratios, and colors in vari- 
ous filters 

— evolutionary and cosmological corrections including 
the effect of attenuation by intergalactic hydrogen. 



Ejection rates of gas and metals by dying stars are cal- 
culated from the individual stellar gas and heavy element 
yields for the respective isochrone metallicities. 

Every star of mass m on an isochrone of age t and 
metallicity Z is attributed a spectrum (&(A, m)) from 
Lejeune et al. 's spectral library, appropriate for its effec- 
tive temperature T e g and surface gravity g. Because our 
model is a 1-zone model without any spatial resolution, 
spectra of all stars along an isochrone of given metallicity 
Z can be summed up, weighted with the IMF $(m), to 
give the integrated isochrone spectrum (S , ; so (A, t, Z)), i.e. 
the spectrum of an SSP of metallicity Z at age t 

/>m u 

S iso {X,t,Z)= S*(\,m,t,Z) ■ $(m) • dm. (1) 

J mi 

Alternatively, the code can also model arbitrary stel- 
lar populations with extended or complex star formation 
histories, such as galaxies. For the evolution of a galaxy 
with a star formation rate extended in time (SFR(t)), 
isochrones are followed as a function of time and multi- 
plied with the star formation rate at the formation time 
t\ of the respective isochrone, SFR(t\). So the integrated 
spectrum of a galaxy (S ga \(\, t, Z)) at time t is calculated 
via Eq. (§) 



Sgal(A, t, Z) 



Siso(A, [t - t{), Z) ■ SFR{t\) ■ dt\. (2) 



To account for the fact that, in general, the stellar 
population of a galaxy is comprised of subpopulations of 
different metallicities, an additional integration over all 
the metallicities can be performed. Here, however, we only 
present results for our SSP models. 

2.5. Filter systems 

Several filter systems are implemented in our evolution- 
ary synthesis code, as e.g. Johnson U B V R I (cf. Lamia 
19820 , J H K L (Bessel & Brett |1988| ), g r (Th uan 



Gunn |l976| ), U+ J+ F+ N+ (Kron p8fj[ Koo p985| ) , HST 
F300W . . F 8UW, Washington C M Tl T2 (Harris & 
Canterna 1979| ), and Stromgren u v b y (cf. Lamia 1982). 
Results for all of them are presented in our electronic ta- 
bles. 

Calibrations for the Washington filters were done with 



colors for stars in Gunn & Stryker's (1983) spectral li- 
brary (Geisler 2000 priv. comm.). The Johnson, Thuan & 
Gunn and Bessel & Brett filter systems were calibrated 
with information from Landolt-Bornstein VI/2b. The 
HST-Filters were calibrated with Kurucz's Vega model 
as described on the HST home page. The calibration 
of Str omgren photometry was done according to Gray 
( |1998|) . 

With the time evolution of SSP spectra as presented 
here, any other filter system within the long wavelength 
range of our spectra from 90 A through 160 /i m is readily 
applicable. Integrating over the spectra with filter trans- 
mission and detector response functions allows to study 
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Table 1. Evolution time t and BDM as a function of 
redshift z with the cosmology Ho = 75, fio = 0.1, Ao = 0.9 
on the left, and with Ho = 65, Qrj = 0.1, Aq = on the 
right hand side, respectively. 



z 


t [Gyr] 


BDM 


2 


t [Gyr] 


BDM 


0.0 


15.12 


36.512 


0.0 


11.92 


31.15 


0.5 


9.932 


42.487 


0.5 


6.957 


42.267 


1.0 


6.720 


44.375 


1.0 


4.535 


44.130 


1.5 


4.632 


45.543 


1.5 


3.115 


45.305 


2.0 


3.222 


46.377 


2.0 


2.189 


46.179 


2.5 


2.236 


47.019 


2.5 


1.542 


46.877 


3.0 


1.522 


47.536 


3.0 


1.065 


47.461 


3.5 


0.988 


47.968 


3.5 


0.702 


47.961 


4.0 


0.579 


48.337 


4.0 


0.417 


48.400 


4.5 


0.258 


48.659 


4.5 


0.188 


48.790 


5.0 


0.0 


48.945 


5.0 


0.0 


49.141 



both the time evolution and the mctallicity dependence of 
any arbitrary filter system. 

2.6. Cosmological models 

As outlined in Sect. [I], our SSP results are not only in- 
tended for use in the interpretation of star cluster data but 
also for combination with any kind of dynamical galaxy 
formation and/or evolution model. For convenient use in 
this second respect, we also calculate evolutionary and 
cosmological corrections. 

Spectra are redshifted to X'(z) = A • (1 + z) and fluxes 
are decreased by F' x (X'(z), z) = F x (X'(z),0)/(l + z). 

For galaxies evolutionary and cosmological corrections 
are conventionally given in terms of magnitude differences, 
which are normalized to the average luminosity of locally 
observed galaxies of the respective type. SSPs, however, do 
not have such an absolute magnitude scale. We therefore 
chose to give their evolutionary and cosmological correc- 
tions in terms of luminosity ratios rather than magnitude 
differences and call them e\ and K\, respectively, not to 
be confused with the e\ and k\ conventionally given for 
galaxies: 

_ L x (z,t(z)) 
L\(z,t ) 
_ L\(z,t Q ) 

Kx ~ L X {o,t y 

with t(z) being the age of a galaxy observed at red- 
shift z and to the age of a galaxy at 2 = 0. Table 
[l] gives the ages t(z) of a stellar population formed 
at redshift Z{ = 5.0 and the Bolometric Distance 
Modulous (BDM(z)) for two specific cosmological mod- 
els (H , n , A ) = (75, 0.1, 0.9) and (65, 0.1, 0) we 
chose to present here. With the large number of spectra 
given for every SSP(Z) in short time intervals any arbi- 
trary cosmological model can easily be applied. 



3. Spectral Evolution of SSPs for 5 different 
metallicities 

We present the time evolution of SSP spectra of 5 metal- 
licities from 4-10 6 yrs (which is the youngest isochrone) up 
to 14 Gyr with a Salpeter and a Scalo IMF, respectively. 
The initial stellar mass is 1.6 • 10 9 M for both IMFs. 

3.1. Time evolution of spectra 
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5000 10000 15000 20000 25000 
MA] 

Fig. 2. Spectrum in terms of flux F\ as a function of wave- 
length A at 3 different times for an SSP of solar metallicity 
and Salpeter IMF. 

In Fig. H we present the spectral evolution of an SSP 
with solar metallicity and Salpeter IMF at 3 different 
ages of 1, 5, and 10 Gyr, respectively. The flux decreases 
with time as stars of lower and lower masses die with the 
strongest effect seen between 1500 and 5000 A. Note the 
small rise of the flux below 1000 A due to very hot stars 
on the white dwarf cooling sequence. 

The 1 Gyr spectrum shows clearly the Balmer, 
Paschen, and Bracket series of hydrogen absorption lines. 
Note the almost unchanging flux at A < 1000 A and the 
strong time evolution in the 4000 A break. 

The spectral evolution of an SSP with a Scalo IMF 
looks very similar to the one with a Salpeter IMF plot- 
ted here, because, at 1 Gyr, most of the high mass stars 
that were less numerous in an SSP with Scalo IMF, have 
already died. A small difference in absolute flux would be 
seen due to the higher number of low mass stars in the 
Scalo IMF, which dominate the flux at higher ages. 

3.2. Metallicity effects on the spectral evolution 

In Figs. H a) and b), we show spectra for SSPs of 3 differ- 
ent metallicities at ages of 1 and 10 Gyr, respectively, all 
with a Salpeter IMF. Enormous differences are seen both 
at young and old ages between SSPs of different metallici- 
ties. As compared to a solar metallicity SSP, the spectrum 
of a low metallicity SSP shows significant differences both 
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3.3. Comparison with earlier models 

We compare our results with those from Bruzual & 
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Fig. 3. Flux F\ as a function of wavelength for SSPs with 
Salpeter IMF and 3 different metallicities at ages of a) 1 
Gyr and b) 10 Gyr. (Z = 0.0004 = 0.02 Z®, Z = 0.008 = 
0.4 Z Q , Z = 0.02 = Zq). 



in continuum slope and absorption line strengths. Over 
the UV and optical wavelength range the continuum spec- 
trum of a low metallicity SSP has stronger flux and the 
maximum at shorter wavelengths as compared to a solar 
metallicity SSP. This is a consequence of stars - on aver- 
age - being brighter and hotter at lower metallicity. Note 
the very strong UV flux at 1000 A < A < 2000 A of the 
1 Gyr very low metallicity (Z = 0.02 Z Q ) SSP. Not only 
the metal absorption lines are weaker - as expected in 
low metallicity SSPs of all ages, but also the Balmer ab- 
sorption lines, most prominent in the 1 Gyr old SSP s, are 
visibly stronger at low metallicities (cf. Kurth et al. 1999 
for a quantitative analysis of absorption features). 

At young ages ~ 1 Gyr, the near infrared (NIR) con- 
tinuum is largely metallicity independent while at older 
ages, e.g. at 10 Gyr as shown in Fig. |3|b, it is visibly af- 
fected by metallicity effects. E.g. in the 1.7 /xm bump in 
the H— band, and at 2.2 /xm in the if— band a 10 Gyr old 
Z = 0.02 Z Q SSP only shows ~ | of the flux of a solar 
metallicity SSP. 



Chariot's ( 1993 ) models (hereafter BC) as obta ined f rom 
the CD-ROM accompanying Leitherer et al. ( 1996 ). In 
Fig. |] we show a comparison of spectra from 1000 A to 
25000 A a), a more detailed view from 5000 A to 10000 
A b) and a view from to 2000 A c) for SSPs with so- 
lar metallicity and Salpeter IMF. Both models show good 
agreement from the UV to ~ 7000 A. Between 7000 A and 
12000 A, however, BC's spectrum shows significantly less 
flux than ours. The same happens again in the H— and 
if— band regions. This probably is due to the TP-AGB 
phase which involves stars with an effective temperature 
around 3000 K. The comparison in the wavelength range 
from 90 A to 2000 A shows a strong difference between 
the models. In this wavelength region very hot stars dom- 
inate the flux. These are stars on the white dwarf cooling 
sequence. We do not know in detail how theses stars are 
handled in BC's models but we suspect that the flux dif- 
ference is due to a different treatment of white dwarfs. 

4. Luminosity Evolution, Colors, and 
Mass-to- Light Ratios 

4.1. Luminosity evolution 

In Fig. ^| we present the time evolution of V— , J—, and 
if— band absolute magnitudes of SSPs with 3 different 
metallicities. Again, the metallicity effects are discussed 
on the example of a Salpeter IMF. 

All over the time evolution, SSPs are brighter at lower 
metallicities in all bands U, . . . , I (Fig. ||a). The differ- 
ence in My between a solar metallicity SSP and one with 
0.02 Zq increases from < 0.5 mag at few 10 s yr to ~ 1 
mag at 5 - 14 Gyr. 

As seen in Fig. ||b, the time evolution of Mj does not 
show any significant metallicity dependence, as could al- 
ready be expected from the spectra in Figs. |^a, b. The 
spectra also indicate what is clearly found for the H— and 
K— band evolution: the trend of luminosity with metallic- 
ity is reversed in the near infrared spectral region as com- 
pared to optical wavelengths with the higher metallicity 
SSPs being brighter in H and K (cf. Fig. ||c), although 
the luminosity differences are not as strong as they are at 
shorter wavelengths. 

4.2. Color evolution 

The (V — I c ) and (B — V) color evolution for SSPs of var- 
ious metallicities with a Salpeter IMF is plotted in Fig. 
||. Note that I c now is the Cousins I for the purpose of 
comparison with data for M31 clusters. On the right hand 
side of Fig. ^ we give the median colors with their 1 a 
dispersions of M31 clusters which have a mean metal licity 
[Fe/H] = -1.2 as determined by Barmby et al. (|2000|) . We 
note that observed integrated GC colors are in good agree- 
ment with our models for the corresponding metallicity at 
ages 12 - 14 Gyr. 
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Fig. 4. Spectra in terms of flux as a function of wave- 
length. Comparison between our model and the model 
from Bruzual & Chariot at 12 Gyr with Salpeter IMF 
and solar metallicity in 3 different wavelength ranges. 
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Fig. 5. Time evolution of absolute magnitudes in Johnson 
V— (a), J- (b), and if -bands (c) for SSPs with 
Salpeter IMF and 3 different metallicities. (Z = 0.0004 = 
0.02 Z Q , Z = 0.008 = 0.4 Z Q , Z = 0.02 = Z Q ). 



During the first ~ 1 Gyr, (V — I c ) colors of all three 
SSPs show a complex behavior with the solar metallicity 
SSP starting out ^0.1 mag bluer and reddening much 
more than the one with 0.02 Z Q . For all ages > 1 — 2 
Gyr, however, SSPs of different metallicities show largely 
different (V"— / c ) colors. E.g., at an age of 12 Gyr, the SSPs 
with Z Q , 0.4 Z Q and 0.02 Z Q have (V - I c ) = 1.3, 1.17, 



and 0.87, respectively. Hence, if e.g. star cluster ages are 
known to be > 2 Gyr, their observed (V — I c ) colors allow 
for largely unambiguous metallicity determinations. 

Note that differences between the two IMFs are very 
small, as already mentioned above, and therefore not 
shown here. In particular, they almost do not affect age 
and/or metallicity determinations at ages > 4 • 10 s yr. 
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Fig. 7. Time evolution of the total stellar mass M in units 
of 10 9 M Q for SSPs with a) Salpeter and b) Scalo IMF and 
three different metallicities. (Z = 0.0004 = 0.02 Z Q , Z = 
0.008 = 0.4 Z Q , Z = 0.02 = Z Q ). 



4.3. Mass-to-light ratios 

In Fig. fj] we show the time evolution of the total mass in 
stars, M, in SSPs with Salpeter and Scalo IMFs and three 
different metallicities. In all cases, SSPs started out with 
an initial stellar mass M = 1.6 • 10 9 M s as mentioned 
above. 

The total stellar mass M decreases as stars lose mass 
by stellar winds and die after ejecting mass either in a PN 
or in a SN event. With a Salpeter and a Scalo IMF, SSPs 
lose ~ 19 % and ~ 23 %, respectively, of their initial mass 
until an age of 12 - 14 Gyr. The bulk of the mass loss 
occurs during the first ~ 2 Gyr in case of a Salpeter IMF 
due to its larger number of high mass stars, while mass 
loss occurs more steadily with a Scalo IMF. 

While for an SSP with Salpeter IMF the metallicity 
dependence of the mass loss is very small, it is slightly 
stronger in the case of a Scalo IMF (23% for Z and 20% 
for 0.02 Z Q ). 

Although the difference in mass loss is small between 
SSPs with Salpeter and Scalo IMFs, the differences in 
mass-to-light ratios M/L\, with A denoting any filter 



band, are much larger because of the luminosity differ- 
ences. 

Fig. U shows the time evolution of V— band mass-to- 
light ratios M/L v of SSPs with both IMFs and three 
metallicities. The mass is the time evolving stellar mass 
of the SSP with mass loss included as depicted in Fig]?]. 

Note the difference of more than a factor of 2 in M/L\ 
at all times and for all metallicities between a Salpeter and 
a Scalo IMF and also the strong metallicity dependence 
in both cases. 

At 10 Gyr, a Salpeter SSP has M/L v - 12, - 9, - 5 
in solar units for Z — Z Q , 0.4 Z Q , and 0.02 Z@, respec- 
tively. For an extensive discussion of M/L— values (and 
their uncertainties) both theoretically derived by various 
groups using different sets of input physics and in compar- 
ison with observationally derived M/L— values obtained 
f or GC s using various methods see Fritze - v. Alvensleben 
( 120001 ). 
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5. Theoretical calibrations of colors in terms of 
metallicity and comparison with GC data 

When spectroscopic abundance determinations are inac- 
cessible or exceedingly time-consuming, like e.g. for GCs 
in external galaxies, broad band colors in combination 
with empirical calibrations with well known abundances 
are used (cf. Couture et al. 1990p . 

Clearly, some colors and color systems (like e.g. the 
Stromgren system) are better suited for this purpose than 
others. It is important to note, however, that in particular 
for optical and NIR colors the relations between color and 
metallicity are more or less age-dependent. This is seen 
in Figsga), b), where we plot theoretical calibrations ob- 
tained from our SSP models at various ages for (B — V) 
(Fig.|a)), (V-I c ) (Fig.gb)) and (V-K) (Fig. go)) versus 
[Fc/H]. We also include observations of GCs in the Milky 
way f rom the GC catalogue compiled by W. E. Harris 
(|1996|) (we only use clusters with E(B — V) < OA ) and 
median colors of GCs in M31 from Barmby et al. ( 2000 ) 
in Fig. ^ which show good agreement with our models 
at 10 — 14 Gyr and low metallicities. In Fig. ^c) observa- 



tions from Brodie & Huchra (1990) are included which also 
show very good agreement with our models. In this low 
metallicity regime empirical calibrations can be used with 
good accuracy. In all three colors it is seen that at young 
ages the theoretical calibrations are significantly different 
from what they are at old ages. Strong changes are seen 
between 0.5 and 1 Gyr and between 1 and 5 Gyr, smaller 
changes at older ages, e.g. between 5 and 14 Gyr. Note 
that theoretical calibrations can reach the highest metal- 
licities for which stellar evolutionary tracks or isochrones 
are available, i.e. 2.5 Zq in our case, while observational 
calibrations are restricted to the metallicity range of Milky 
Way GCs, i.e. to Z < 0.3 Z & . 

Over the metallicity range of the GC sample —2 < 
[Fe/H] < —0.5 the agreement between our theoretical cal- 
ibrations at 12 Gyr and the empirical ones from Couture 
et al. p90| ) is very good, i.e. better than 0.05 mag both 
in (B-V) and (V-I c ). 

Note, however, that extrapolations from the approx- 
imately linear observational relation at [Fe/H] < —0.5 
to higher metallicities would produce grossly misleading 
results since towards higher metallicities the theoretical 
relations become significantly steeper. By extrapolating 
the empirical linear relation, a 10 — 14 Gyr old cluster 
with e.g. (B — V) = 1.1 would be ascribed a metallicity 
Z cx tia,p ~ 4.4 Zq while our model shows it to only have 
Z w 2.5 Zq. 

For (V — I c ), the deviation from the linear relation 
valid in the low metallicity region is even stronger. E.g. 
an old cluster with (V — I c ) — 1.2 would be ascribed 
•^extrap ~ 3.0 Zq while our models indicate solar metal- 
licity, i.e. three times lower! We stress that in particular 
for investigations of intermediate age cluster populations 
or of clusters with metallicities that might exceed those of 
the Milky Way or M31 GCs (on which the empirical cal- 
ibrations are based) it is very dangerous to use empirical 
calibrations. Our models clearly show the age dependence 
as well as the strong discrepancies from the empirical lin- 
ear color - metallicity relations for [Fe/H] > —0.5. 

In Figs. ^2 a), b) we present analogues to Figs. ^|a), b), 
c) for the Stromgren color indices mi — (v — b) — (b — y) 
and ci = (u — v) — (v — b). For the particularly metal- 
sensitive index mi it is seen that it allows to well separate 
clusters in terms of metallicities and ages. It allows for a 
better metallicity separation on the basis of colors than 
the Johnson filter system does. But it has a smaller ab- 
solute difference, which makes it fragile to observational 
uncertainties. 

The age-sensitive Stromgren color index c\ shows a 
completely different behavior. While it is, indeed, rela- 
tively constant at given age for all [Fe/H] < —0.5, it starts 
changing significantly towards higher metallicities, where 
it allows for metallicity separations nearly as good as m\. 
Due to its changing behavior, it is important to well know 
the metallicity of the cluster, which should be < —0.5 in 
[Fc/H]. Above this metallicity, even mi shows a better age 
separation than c\. 
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While empirical calibrations like those obtained from 
Milky Way or M31 GC analyses are only applicable for 
objects of comparatively old ages, our theoretical calibra- 
tions can also be used for intermediate age and young clus- 
ters, as long as emission lines from the active SF phase do 
no longer affect the broad band colors. From our exten- 
sive electronic tables, theoretical calibrations analogous 
to the ones depicted in Figs [)]a), b), c) can be extracted 
for any color of our filter systems. They provide a useful 
tool for the interpretation of young star clusters forming 
in large numbers in interacting, merging, and starburst- 
ing gas-rich galaxies, for young and intermediate age GCs 
in old merger remnants and kinematically young ellipti- 
cal and SO galaxies. The extension to mctallicities in the 
range 0.5 Zq . . . 2.5 Zq, that is only possible with theo- 
retical calibrations, is of particular importance since spec- 
troscopy of the first handful of bright young star clusters 
recently formed in spiral galaxy mergers has proven their 
relatively high metallicity - co mpared to Galacti c GC s 
(e.g. Schweizer & Seitzer |1993| , Whitmore et al. |l999p . 
Since these 2 nd generation GCs are formed from gas pre- 
enriched in the interacting and/or starbursting galaxies, 
metallicities > 0.5 Zq were expected (Fritze - v. A. & 
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Table 2. Evolutionary and cosmological corrections e\ 
and K\, respectively, as a function of redshift z from z = 
to z — 5 in the cosmology Hq = 65, f2o = 0.10 for 
all Johnson bands for an SSP with solar metallicity and 
Salpeter IMF. Full table with all wavelengths bands and 
for all metallicities and two different cosmologies is given 
in the electronic version, at CDS, and on our web page. 



today's age to = t(z = 0) via 

m x {z) = -2.5log(L x (0,t )-e\(z)-K X (z))+BDM(z) 



6. Cosmological Evolution 

In this section, we present evolutionary and cosmological 
corrections for direct use of our SSP results in combination 
with galaxy formation or cosmological structure formation 
scenarios. We chose two different cosmological models as 
characterised by (Ho, Hq, A ) = (65, 0.1, 0.0) and 
(75, 0.1, 0.9), respectively, and an assume an onset of SF 
at redshift Zf = 5. The manifold of cosmological models is 
constrained by the ages of Milky Way GCs (12 - 14 Gyr). 

In case of galaxies, the luminosities of model galaxies 
at today's age of 12 - 14 Gyr are normalized to their ob- 
served average local luminosities, i.e. to the average local 
(Mb) of the respective galaxy type, and their evolution- 
ary and cosmological corrections are then given in terms of 
magnitude differences. Such a normalization is not mean- 
ingful for SSPs if understood as galaxy constituents or 
subpopulations. Hence, as explained in Sect. |2.6| , we chose 
to present their evolutionary and cosmological corrections 
(e x , k\) in terms of luminosity ratios in selected filter 
bands A. 

In Table || we present these corrections as a function 
of redshift z for the Johnson U . . . K bands. They are 
given for the two cosmological models chosen above. Only 
the first few columns and lines of this table are shown for 
one cosmology in the printed version for orientation. Full 
tables are available in the electronic version, at CDS, and 
our web site. 

The apparent magnitude m\(z) in some filter A of 
an SSP of given mass and metallicity at redshift z, 
m\(z) = — 2.5 log L\(z,t(z)) is obtained from the 
absolute luminosity L x (0,to) of an SSP of given mass by 



7. Conclusions 

A new set of evolutionary synthesis models are presented 
for single burst single metallicity stellar populations cov- 
ering metallicities 0.02 < Z/Zq < 2.5 and ages 4- 10 6 yr < 
t < 16 Gyr. They are based on the most recent isochrones 
from the Padova group (Nov. 1999) that extend earlier 
models by the inclusion of the thermal pulsing AGB phase 
for stars in the mass range 2 Mq < m < 7 Mq in ac- 
cordance with the fuel consumption theorem. We show 
that with respect to earlier models, inclusion of the TP- 
AGB phase leads to significant changes in the (V — /) 
and (V — K) colors of SSPs in the age range from 10 s 
to > 10 9 yr, in particular at metallicities Z > 0.5 Zq. 
The ages derived from observed (V — I) colors of young 
star clusters in this age range decrease by a factor of ~ 2 
when the TP- AGB phase is includ ed. Using m odel atmo- 
sphere spectra from Lejeune et al. ( 1997 , 1998 ), we calcu- 
late the spectral evolution of single burst populations of 
various metallicities covering the wavelength range from 
90 A through 160 /ini. Expanding dust shells around 
AGB stars, however, are not accounted for by our models. 
Therefore, they should not be used for the interpretation 
of IR data (i.e. beyond K) of intermediate age star clus- 
ters. Age and metallicity effects of the spectral and photo- 
metric evolution of our SSP models are discussed in detail. 
Isochrone spectra are convolved with filter response func- 
tions to describe the time evolution of luminosities and col- 
ors in Johnson, Thuan & Gunn, Koo, HST, Washington, 
and Stromgren filters. 

Theoretical calibrations are provided for a number of 
colors in terms of metallicity and compared with empiri- 
cal calibrations and data from GCs. It is found that the 
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empirical calibrations are only suitable for the age and 
metallicity range of the clusters they are derived from. 
Our models do not have such restrictions and cover the 
whole metallicity - and age - range. They allow to study 
color - metallicity calibrations in their time evolution and 
up to metallicities > Zq, both of which is important for 
the interpretation of young and intermediate age star clus- 
ter systems. 

Model results are not only intended for use in the inter- 
pretation of star clusters over the full age range from few 
10 — 1.4-10 10 yr but also for combination with any kind of 
dynamical galaxy formation and/or evolution model that 
contains a star formation criterion, i.e. some prescription 
how to transform gas into stars. Moreover, the evolution 
of these single burst single metallicity stellar populations 
is readily folded with any kind of star formation - and 
eventually chemical enrichment - history to describe the 
evolutionary spectral synthesis of composite stellar pop- 
ulations like galaxies of any type with continuous or dis- 
continuous SF. 

For these latter purposes we also present the time evo- 
lution of mass-to-light ratios (including ejection rates for 
gas and metals) for two different IMFs as well as cosmo- 
logical and evolutionary corrections for all the filters as 
a function of redshift z for < z < 5 and two different 
cosmologies. 

Extensive data files are provided in the elec- 
tronic version, at CDS, and at http://www.uni- 
sw.gw jg.dc/^galev/ . 
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Fig. 10. a) Stromgrcn mi and b) c\ color indices vs. 
metallicity [Fe/H] for SSPs with Salpeter IMF at various 
ages. 0.5 Gyr (open circles), 1 Gyr (filled squares), 5 Gyr 
(stars), 10 Gyr (open squares), and 14 Gyr (filled circles). 



Fig. 9. a) (B — V), b) (V - J c ) (J c is Cousins I) and c) 
(V-K) color vs. metallicity [Fe/H] for SSPs with Salpeter 
IMF at various ages. 0.5 Gyr (open circles), 1 Gyr (filled 
squares), 5 Gyr (stars), 10 Gyr (open squares), and 14 
Gyr (filled circles). Crosses in a) and b) are Milky Way 
GC data from Harris' catalogue. In c) the crosses are M31 
GC data from Brodie & Huchra. Also included in all three 
figures are the mean M31 GC data and their la ranges 
from Barmby ct al. . 



